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Abstract 

We study the charge carriers transport in an AA-stacked bilayer graphene modulated by a 
lateral one-dimensional multibarrier structure. We investigate the band structures of our system, 
that is made up of two shifted Dirac cones, for finite and zero gap. We use the boundary conditions 
to explicitly determine the transmission probability of each individual cone (r = ±1) for single, 
double and finite periodic barrier structure. We find that the Klein tunneling is only possible when 
the band structure is gapless and can occur at normal incidence as a result of the Dirac nature 
of the quasiparticles. We observe that the band structure of the barriers can have more than one 
Dirac points for finite periodic barrier. The resonance peaks appear in the transmission probability, 
which correspond to the positions of new cones index like associated with r = ±1. Two conductance 
channels through different cones (r = Tl) are found where the total conductance has been studied 
and compared to the cases of single layer and AB-stacked bilayer graphene. 
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1 Introduction 


Graphene, single layer of pure carbon atoms in a honeycomb lattice, has a gapless linear electronic 
spectrum near to the two valleys of the first Brillouin zone [1-3]. That yields electrons in graphene are 
described by a massless two dimensional relativistic Dirac equation and the corresponding electrons 
have a chiral nature [4,5]. This leads to many unusual electronic properties and potential applications 
[2,3]. Another interesting property is the Klein tunneling [ 6 ], describing that the Dirac fermions can 
be transmitted with probability one through a classically forbidden region. Graphene can not only 
exist in the free state, but two or more layers can stack above each other to form what is called 
few-layer graphene, as example bilayer graphene is resulting from a composition of two stacked sheets. 
Bilayer graphene has surged as another attractive two-dimensional carbon material. In addition, it is 
demonstrated that the bilayer graphene has a new unusual physical properties and the spectrum is 
different from that of single layer graphene. There are two dominant ways in which the two layers can 
be stacked: the AB-stacked and AA-stacked bilayer graphene. 

AB-stacked bilayer graphene or Bernal form is the basis of the graphite from [7] and is usually 
derived. In which only one atom in the lower layer lies directly below an other atom in the upper layer 
and the other two atoms over the center of the hexagon in the other layer. AB-stacked bilayer graphene 
has a gapless quadratic dispersion relation, two conduction bands and two valance bands, each pair is 
separated by an interlayer coupling energy of order 71 = 0.4 eV. Recent experiments showed that the 
AA-stacked bilayer graphene could also exist and is a new stable stacking of graphene [8-13], where 
the A sublattice of the top layer is stacked directly above the same sublattice of the bottom layer. 
In AA-stacked bilayer graphene, the energy bands are just the double copies of single layer graphene 
bands shifted up and down by the interlayer coupling 7 = 0.2 eV. Additionally, it has band structures 
that differ from those corresponding to single layer and AB-stacked bilayer graphene. 

Recently, many works have been reported in order to investigate the multibarrier structures in 
single layer graphene [14-21]. However, much less experimental and theoretical works have been done 
on AB-stacked bilayer graphene [ 21 , 22 ]. Motivated by different developments on the subject and in 
particular [14,23,24], we investigate the energy bands for finite and zero gaps in AA-stacked bilayer 
graphene. Our theoretical model is based on the well established tight binding Hamiltonian [23-26] 
and we adopt the parametrization of the relevant intralayer and interlayer couplings. We note that 
the electronic band structure can be modified by the application of a periodic potential and the gaps. 
We explicitly calculate the transmission probability of each individual cone (r = ±1) for single, double 
and finite periodic barrier structure. We show that the Klein tunneling is only possible when the band 
structure is gapless and can occur at normal incidence. This makes difference with respect to the case 
of AB-stacked bilayer graphene where no Klein tunneling (T = 1) is expected [21,27]. We show that 
the existence of the intracone transition allows only two transmission channels in AA-stacked bilayer 
graphene resulting in a total conductance, while the transmission in the intercone is zero. In contrast, 
for AB-stacked bilayer graphene [27] we have intracone and intercone transitions, which are possible 
between all four bands. In addition, the total conductance of AA-stacked bilayer graphene is different 
with that of the single layer graphene [14]. 

The rest of the paper is organized as follows. In section 2, we consider the multibarrier structure 
for the AA-stacked bilayer graphene. We formulate our model by setting the Hamiltonian system 
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and calculating the associated energy bands in each potential region. We obtain the spinor solution 
corresponding to each regions. Using the transfer matrix at boundaries together with the incident, 
transmitted and reflected currents we end up with two transmission probabilities. In section 3, we 
numerically present our results of the transmission probability of each individual cone for single, 
double and flnite periodic barrier structure. In section 4, based on the obtained results we study two 
conductance channels through different cones and as well as the corresponding total conductance. We 
conclude our work in the flnal section. 

2 Energy bands 

We consider electron tunneling through a lateral one-dimensional multibarrier in AA-stacked bilayer 
graphene as shown in Figure 1(c), consisting of + 2 regions indexed by r. The parameters for 
the regions represented by an odd index r = 2 n_ + 1 , are the potential barrier heights V2n-+i^ the 
widths d and the gaps A 2 n_+i- However, for an even index, we have V 2 n^ — — 0, where 

n± = {O, 1, • • • , with the widths a. 



Figure 1: (a): Illustration of lattice structure of AA-stacked bilayer graphene, which is consisting of two 
coupled single layer graphene, (b): unit cell of AA-stacked consists of four atonns each pair is separated by 
an interlayer coupling energy of order 7 = 0.2 eV. (c): nnultibarrier structure consisting of A -h 2 regions. 


The charge carriers in the AA-stacked bilayer graphene are described, in each region r, by the following 
four-band Hamiltonian [23] 

^ Vr + Ar VpiPx-iPy) 7 0 ^ 

^ _ '^f{,Px ^Py) V Aj. 0 7 
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The eigenstates of are the four component spinors j , where and 

^ 5 (s') envelope functions associated with the probability amplitudes of the wave functions 

on the A{A') and B{B') sublattices of the lower (upper) layer. In (1), px^y = —ihVx^y is the in-plane 
momentum relative to the Dirac point, — lO^m/s is the Fermi velocity for electrons in single layer 
graphene, Vr is the potential heights and is the gaps. In Figure 1(a), we have the bilayer graphene 
consisting of two layers of graphene having the structure AA-stacked, where the A sublattice of the 
top layer is stacked directly above the same sublattice of the bottom layer by the interlayer coupling 
7 = 0.2 eV. The unit cell is consisting of four atoms labeled A, B on the lower layer and A\ B' on 
the upper layer, see Figure 1(b). 

We notice that our system is infinite along the ^-direction since [Hr^Py] = 0 and then we can 
decompose our spinor into 

^^{x,y) = e^^yyr{x,ky). (2) 

From the eigenvalue equation we end up with four linear differential equations of the 

from 


-irj(d^ + ky)ip^s + = {E-Vr- Arm 

(3a) 

-iv{d. - kyWA + Vb' = {E-Vr + ArWB 

(3b) 

+ kyWe' +rA = iE-Vr- ArWA' 

(3c) 

-irjidx - ky)m +Vb = {E -Vr + Ar)^ 

(3d) 


where we have introduced the length scale p = « 3.29 nm. From now on we switch to dimensionless 

quantities by measuring all energy terms in units of the interlayer coupling 7 , namely E — 

Vr —^ and We combine the above equations to eliminate the unknown ones at time 

to end up with the following differential equation 

[dl +{Kf]'fpB{x,ky) = 0 (4) 


where the wave vector along the a;-direction reads as 

K = ^-kl + y-^{{E-Vr-Tf-^i) ( 5 ) 

and T is the cone index, with r — —1 {r — +1) for the lower (upper) cone. Then we have V2n-+i and 
^2n_+i equal to Vi and Ai, respectively. For the regions where Vi = Ai = 0, we find the wave vector 

k^ = (6) 

and the corresponding energy band is 

E^^’- = T + s^^r,^{{klY + {kyY). ( 7 ) 

However generally, for any region we can deduce the energy from previous analysis as 

E^>^ = V + r + {{klY + {kyf) + Af (8) 

where si is the chirality index of a quasiparticle in the barrier regions (Vr = Vi, A^ = Ai) and in 
outside (Vr = A^ = 0) is sq, with si = +1 (or sq = +1) and si = —1 (or sq = —1) are the electron-like 
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and hole-like band indexes. It is important to note that in ( 8 ), for 7 —^ 0, the energy bands will be 
reduced to 

E^ = Vi + sr^hvF{{klY + {kyY) + A? (9) 

which is a degenerate energy and corresponding to two monolayers. It is clear that this degeneracy 
can be lifted by taking into account the coupling between the two monolayers, i.e. 7^0 [28]. 



-0.4 -0.2 0.0 0.2 0.4 
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Figure 2 : The energy bands as a function of the nnonnentunn ky, where the blue line (green line) corresponds 
to the lower cone (upper cone), (a): for Vi = Ai = 0. (b): for Ai = 0 and Vi = (7, I.67, 2.47) 
correspond to (dashed line, solid line, dotted line), (c): for Vi = 2.47 and Ai = (O.27, O.47, O.87) 
correspond to (dashed line, solid line, dotted line). 


In Figure 2, we plot the energy bands as a function of the momentum ky for the lower cone r = — 1 
(blue line) and the upper cone r = 1 (green line). In Figures 2(a) and 2(b) the energy bands are plotted 
for Ai = 0 while in Figure 2(c) for Ai ^ 0. It is evident that in Figure 2 the energy bands are just the 
double copies of monolayer graphene bands, with one Dirac point is shifted to wards positive (r = +1) 
and the other towards negative (r = —1) energies. For zero gap the spectrum is linear (Figures 2(a) 
and 2(b)), however for a finite gap the spectrum is parabolic (Figures 2(c)). In Figure 2(b), we plot 
the energy bands for three different values of the potential heights Vi = ( 7 , I. 67 , 2 . 47 ). From this, 
we can observe that when we increase the potential heights Vi, the energy bands increase upwards. To 
see the effect of the gap, we plot the energy bands for three different values of Ai = (O. 27 , O.dy, O. 87 ) 
and for Vi = 2 . 47 , the results are shown in Figure 2(c). It is clearly seen that a gap appears at each 
of the two Dirac points (Vi + t — Ai < E <Vi + r + Ai). When we increases Ai, the width of these 
gaps increases as well. 

Figure 3 represents the band structure of AA-stacked bilayer graphene for zero and finite gap. 
For r = 0, which corresponds to the incident region, the electron states can be subdivided into two 
regimes. The first one is where both r = ±1 electrons are electron-like, the second one is where 
r = +1 electrons are hole-like while r = — 1 electrons are electron-like. However, for the transmission 
region (r = iV -h 1), we have two regimes as we shall show in the incident region. Note that we have 
four transmissions that are related to the band structures on the incident and transmission region. In 
addition, we notice that for the intracone transitions (i.e. r —^ r processes), there exist four cases 
for transitions across the barrier: 

• Electron in regime I —^ Electron in regime I (1) 

• Electron in regime I —Electron in regime II (2) 
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• Electron in regime II —^ Electron in regime I (3) 

• Electron in regime II —^ Electron in regime II (4) 


However, all intercone transitions (i.e. r ^ —r processes) are strictly forbidden due to the orthog¬ 
onality of electron wave functions with a different cone index [24]. As already mentioned above, the 
transitions depend on the incident and transmission regions, so they depend anywhere in the presence 
or absence of the gap. For each of the four cases for transitions, the transmission can be calculated 
using the same method. 


(a) 




Figure 3: Schematic representation of the band structures, (a): for the barriers structure with zero gap. 
(b): for the barriers structure with finite gap. 


3 Transmission probability 

Next we shall show in detail the calculation of the transmission probability of electrons across the 
periodic barrier structure in our AA-stacked bilayer graphene system. The wave function solution of 
the system in each region r can then be written in terms of a matrix form as 


= Gr- Mr- Ar 


where different matrices are given by 


Gr = 


such that reads as 

and the phase is defined by 


-^^5+ £'T 

Sr'^Jr Sj.TJr 


„iklx 


Mr = 


0 e 


0 

—iklx / ’ 


Jr E-Vl-T-Ai 

(l)^ = arctan(/cy/fcj^). 



( 10 ) 


( 11 ) 


( 12 ) 

( 13 ) 


We are interested in the normalization coefficients, namely the components of A^, on the both sides 
of the multibarrier structure. In other words, for the incident region (r = 0) and transmission region 
{r = N + 1 ), respectively, we have 
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^5 = (l,r"r, ^]v+i = (i",Of (14) 

where and V are the reflection and transmission coefficients of each cone (r = il), respectively. 
We use two different matrix notation written as 


G 2 n+ • = Gq * Mq, G 2 n _+1 * ^ 2 n_+l = • Mi 


(15) 


where we have Vi = Ai = 0 in r = 2 n+ and Vi ^ 0 and Ai ^ 0 for r = 2 n_ + 1 , with n± = 
{O, 1, • • • , Using the boundary conditions and the transfer matrix method, we can connect 

Aq with ^]v+i through the matrix ( 


C — ^ + ^)] 'Gq^ ’Gi ’Ml [j{d + a)] ’M^ ^ [{j + l)d + ja] ’G^ ^ - Go ’^o [(j + + ja] . ( 16 ) 

j=o 


Consequently, we have two channels for the transmission probability in each individual cone, which 
are given by 



(17) 


where ("ii is a element of the matrix 



Figure 4: Density plot of transmission probability as a function of the incident angle 00 and its energy E, 
for Vi = 37 , a = 2(i = 30 nm and Ai = O 7 . (a)/(b) for single barrier with r = —1 /t = 1. (c)/(d) for 
double barrier with r = —l/r = 1. 


To allow for a suitable interpretation of our main results, we compute numerically the transmission 
probabilities for r = 1 and r = — 1. In what follow we start by studying the transmission probability 
for single, double barrier and finite periodic barrier structures. We first consider the case of single 
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and double barrier geometry with zero gap. In Figure 4 we show the density plot of the transmission 
probability as a function of the incident angle and its energy. The different colors from blue to orange 
correspond to different values of the transmission from 0 to 1. It is important to note that in the case 
of AA-stacked bilayer graphene the band structure is composed of two Dirac cones shifted by r = — 1 
(Figures 4(a) and 4(c)) and r = 1 (Figures 4(b) and 4(d)). In addition, the transmission probabilities, 
for T — —1 and r = 1, have the same form as that in the case of single layer graphene [14]. In 
both cases of single layer and AA-stacked bilayer graphene, around the point 00 = 0, we have perfect 
transmission with a manifestation of the Klein tunneling effect [6]. Until now we notice that there 
is some similarity between single layer and AA-stacked bilayer graphene. However, for AA-stacked 
bilayer graphene both the electrons and holes have the same chirality index 5 = ±1 while in the case 
of single layer graphene the electron has always -hi and —1 for the hole. For the transmission of the 
lower layer (Figures 4(a) and 4(c)) the Dirac points correspond to F’ = Vi — 1 and for the upper 
layer (Figures 4(b) and 4(d)) correspond to F = Vi -h 1. The results for a symmetrical double barrier 
structure of width d and interbarrier separation a, are shown in Figures 4(c) and 4(d). Compared to 
the results seen above for single barrier, we observe the appearance of peaks in transmission probability. 





Figure 5: Density plot of transmission probability as a function of the incident angle 0o and its energy F, 
for Vi = 37 , a = 2d = 30 nm and Ai = O.dy. (a)/(b) for single barrier with r = —l/r = 1. (c)/(d) for 
double barrier with r = —1 /t = 1. 

Now let us see what will happen if we introduce a gap in the band structures. To do this, we extend 
the results presented in Figure 4 to the case Ai = O .47 to get Figure 5. It is clearly shown that, in 
contrast to the single barrier case (Figures 5(a) and 5(b)), there are full transmission resonances, in 
the case of double barrier (Figures 5(c) and 5(d)), inside the gap (Vi + r — Ai < E < Vi + r + Ai). 
Note that, the transmission resonances are resulting from the available states in the well between the 
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barriers. As already mentioned above, for a gapless graphene we have perfect transmission that is 
a manifestation of Klein tunneling effect. In addition, it may be noted that the opening gap in the 
barrier region suppresses this effect. 



Figure 6: Density plot of transmission probability as a function of the incident angle 00 and its energy E, 
for finite periodic barrier structure that contain 25 regions, with Ki = 3 and a = 2(i = 30 nm. (a)/(c) for 
r = —1 with Ai = 0 /Ai = O.dy. (b)/(d) for r = 1 with Ai = 0/Ai = O.dy. 

Next we consider the finite periodic barrier structure that contain 25 regions. In Figure 6 we show 
the transmission probabilities, associated with r = ±1, as a function of the incident angle 0o and its 
energy F’, for Vi = 3 and a = 2(i = 30 nm. Figures 6 (a)/(c) for Ai = 0/Ai = O.dy and r = —1 while 
Figures 6 (b)/(d) for Ai = 0/Ai = O.dy and r = +1. It is evident that the band structure of such 
barrier can have more than one Dirac points located at F = Vi — 1 for the lower cone (Figure 6(a)) and 
at F = Vi +1 for the upper cone (Figure 6(b)). It is important to note that the location of Dirac point 
does not depend on the number of barriers. In fact, the position of the Dirac point is the same for 
single, double and the periodic barrier structures. However, increasing the number of barrier leads to 
the occurrence of resonance peaks in the transmission probability. These resonance peaks correspond 
to the positions of new cone index like associated with r = ±1 (Figures 6(a) and 6(b)) in the finite 
periodic barrier structures. We should emphasize that for r —^ 0, our system will be reduced to the 
case of single layer graphene [14,29]. From Figure 6(c) and 6(d), we can clearly see that the presence 
of the gap leads to suppression of the Klein tunneling effect [6]. In addition, it is important to note 
that the number of peaks presented inside the transmission gap (Vi + r — Ai < F < Ki + r + Ai), in 
the case of double barrier structure (Figure 5(c) and 5(d)), is the same as that corresponding to the 
transmission trough finite periodic barrier structure. 



































4 Conductance 


Basing on the obtained results above regarding the transmission probabilities of our system, we will 
see how the two conductance (G^) of each cone channel (r = ±1) will behave. For this purpose, we 
evaluate by using the Landauer-Biittiker [30] formula 

G" = Go / T" {E, 00) cos 00 #o (18) 

Jo 


where the unit conductance 

Go = N Lykpe^ / hiT 


(19) 


with the factor = 4 is due to the spin and valley degeneracy, Ly is the width of the sample in the 
^-direction and 

kF = ^kl + {kiy^sor^-\E-T). ( 20 ) 

The total conductance Gt is defined as the sum of the conductance channels in each individual cone 
(G^) such as [24] 

Gt = - (G+ + G ) (21) 

and the factor 1/2 is required because the total conductance is contribution of the two cones. Whereas, 
in the case of AB-stacked bilayer graphene, there are four transmissions channels, then the total con¬ 
ductance is the sum of four conductances channels [27]. In the forthcoming analysis, we evaluate 
numerically the conductance in AA-stacked bilayer graphene. 
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Figure 7: Conductance in each individual cone as a function of the energy for a = 2d = 30 nm and 
Vi = 3, with single (green line) and double (red line) barrier structure and also periodic barrier contain 25 
regions (blue line). (a)/(c) for r = —1 and Ai = 0/Ai = O. 47 . (b)/(d) for r = 1 and Ai = 0/Ai = O .47 
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We present in Figure 7 the conductance in each individual cone, given in (18), for a single (green 
line), double barriers (red line) and also for a finite periodic barrier contain 25 regions (blue line), as a 
function of the energy F’, for Vi = 3 and a = 2(i = 30 nm. In Figures 7(a) and 7(b) the conductances 
are plotted for Ai = 0 while in Figure 7(c) and 7(d) for Ai = 0.4. It is important to note that, in 
the case of single, double and finite periodic barriers, we have the same position of the Dirac point 
for each individual cone (F^ = Vi — 1 in Figure 7(a) and F’ = Vi + 1 in Figure 7 (b)). We conclude 
that the position of the Dirac point is the same whatever the number of the barriers. In addition, 
for the case of finite periodic barrier we can have more than one Dirac point at the same position. 
Moreover, when we increase the number of barrier the conductance decrease and new peaks appear. 
To see how the gap will affect the conductance channel in each cone (r = ±1), we extend the results 
presented in Figures 7(a) and 7(b) to the case Ai = O.dy to get Figures 7(c) and 7(d). For the case 
of single barrier (green line), the conductance is zero and there are no resonances in the regime of 
energy (Vi + t — Ai < E < Vi + r + Ai). However, for the double barrier structure (red line) new 
peaks of resonances appear in the above mentioned regime of energy. These peaks can be attributed 
to the bound electron states in the well region between the barriers. The number of these peaks, in 
the regime of energy (Vi+r —Ai < E < Vi + t +Ai), is the same as that of the finite periodic barrier. 




E/y 


Ely 


Figure 8: Total conductance as a function of energy, for Vi = 3 and a = 2(i = 30 nm, with single (green 
line) and double (red line) barrier structure and also finite periodic barrier contain 25 regions (blue line), 
(a) for Ai = 0 and (b) for Ai = O.dy 


Now, we plot the total conductance (21) as a function of energy for single (green line), double 
(red line) barrier structure and for finite periodic barrier contain 25 regions (blue line). Indeed, 
the existence of the intracone transition allows only two transmission channels in AA-stacked bilayer 
graphene [24] resulting in a total conductance, while the transmission in the intercone is zero. In 
contrast, for AB-stacked bilayer graphene [27] we have intracone and intercone transitions, which are 
possible between all four bands. In addition, the total conductance of AA-stacked bilayer graphene is 
different with that of the single layer graphene [14]. 
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5 Conclusion 


We have studied the electronic transport of electrons through single, double and finite periodic barrier 
on the AA-staked bilayer graphene, as well as analyzed the transmission probability and corresponding 
conductance we have started by formulating our Hamiltonian model that describes the system under 
consideration and getting the associated energy bands. The obtained bands are composed of two Dirac 
cones shifted up and down by the interlayer coupling. 

Subsequently, we have calculated the transmission probabilities and presented numerically the 
results for each individual cone (r = ±1). For AA-stacked bilayer graphene the transmission in the 
intercone is zero, in contrast, for AB-stacked bilayer graphene where we have intracone and intercone 
transitions. By increasing the number of barriers an occurrence of resonance peaks in the transmission 
probabilities. These resonance peaks correspond to the positions of new cones index like associated 
with T = ±1 in the finite periodic barrier structures. However, the location of Dirac point {E = Vi +r) 
does not depend on the number of barriers. For a gapless graphene we have perfect transmission that is 
a manifestation of Klein tunneling. However, the opening gap in the barrier region leads to suppression 
of the Klein tunneling effect for near-normal incidence. 

Basing on the obtained results for the transmission probabilities we have found that there are two 
conductance channels through different barriers cases. In the case of finite periodic barrier we can 
have more then one Dirac point at the same position. Also, we have investigated the effect of the gap 
on the conductance channel in each cone. It was noticed that in the well region between the barriers 
the appearance of new peaks of resonances, which can be attributed to the bound electron states. 
Finally the total conductance studied as an average of the two conductance channels in each cone and 
we have concluded that this later for AA-stacked bilayer graphene is different with that corresponding 
to single layer graphene. 
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